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Abstract. Magnetic helicity fluxes in turbulently driven dynamos are studied to demonstrate 
their ability to alleviate catastrophic quenching. A one- dimensional mean-field formalism is used 
to achieve magnetic Reynolds numbers of the order of 10^ . We study both diffusive magnetic 
helicity fluxes through the mid-plane as well as those resulting from the recently proposed 
alternate dynamic quenching formalism. By adding shear we make a parameter scan for the 
critical values of the shear and forcing parameters for which dynamo action occurs. For this 
afi dynamo we find that the preferred mode is antisymmetric about the mid-plane. This is also 
verified in 3-D direct numerical simulations. 
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1. Introduction 

The magnetic field of the Sun and other astrophysical objects, hke galaxies, show 
field strengths that are close to equipartition and length scales that are much larger 
than that of the underlying turbulent eddies. Their magnetic field is assumed to be 
generated by a turbulent dynamo. Heat is transformed into kinetic energy, which then 
generates magnetic energy, which reaches values close to the kinetic energy, i.e. they are 
in equipartition. The central question in dynamo theory is under which circumstances 
strong large-scale magnetic fields occur and what the mechanisms behind it are. 

During the dynamo process, large- and small-scale magnetic helicities of opposite signs 
are created. The pre sence of small-scale helicity works against the kinetic a-efFect , whic h 
drives the dynamo (jPouauet et al. Il976l: lBrandenbur3l200ll; [Field fc Blackma"3 12002| ). 
As a consequence, the dynamo saturates on resistive timescales (in the case of a periodic 
domain) and to magnetic field strengths well below equipartition (in a closed domain). 
This behavior becomes more pronounced with increasing magnetic Reynolds number 
Rcm, such that the saturation magne tic energy of the large-scale field decreases with Re^ 



( Brandenburg Sz Subrama,nian 20051). for which it is cal led catastrophic. Such concerns 
were first pointed out by ^ Vainshtein fc Cattaneo (1992). The quenching is particularly 
troublesome for astrophysical objects, since for the Sun ReM = 10^ and galaxies ReM = 
10i«. 



2. Magnetic helicity fluxes 

The first part of this work addresses if fiuxes of small-scale magnetic helicity in an 
dynamo can alleviate the catastrophic quenching. We want to reach as high mag- 
netic Reynolds numbers as possible. C onsequently we consider the mean-field formalism 
( MoffattifigSOt iKrause fc Radleii Il980l) in one dimension, where a field B is split into a 
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mean part B and a fluctuating part b. In mean-field theory the induction equation reads 

dtB ^TjV^B + V X {U xB + £), (2.1) 

with the mean magnetic field B, the mean velocity field U, the magnetic diffusivity t], 
and the electromotive force £ — u x b, where u = U — U and b — B~B are fluctuations. 
A common approximation for £, which relates small-scale with the large-scale fields, is 



£ = aB ?;tV x B, 



(2.2) 



where ?7t = Wrms/(3fcf) is the turbulent magnetic diffusivity in terms of the rms velocity 
Urms and the wavenumber fcf of the energy-carrying eddies, and a = uk + Q-m is the 
sum of kinetic and magnetic a, respectively. The kinetic a is the forcing term, i.e. the 
energy input to the system. In this model ax vanishes at the mid-plane and grows 
approximately linearly with height until it rapidly falls off to at the boundary. The 
magnetic a can be approximated by the magnetic helicity in the fiuctuating fields: am ~ 
h[ X (/iopof?tfcf /-Bgq), where /io is the vacuum permeability, pa is the mean density, B, 



cq 



(MoPo) 



1/2, 



is the equipartition field strength and h{ = a • b the magnetic helicity in 



the large-scale fields. 

The advantage of this approach is that we can use the time evolut ion equation for the 
magn etic helicity to obtain the evolution equation for the magnetic a ([Brandenburg et. al. 
200i) 



da 



M 



dt 



£ B 



Rcm/ 



(2J) 



where JFq is the magnetic helicity flux term. To distinguish this from the algebraic 
quenching ( Vainshtein &: Cattaneol 1992 ) it is called dynamical a-quenching. 

For the flux term on the RHS of equation (|2.3p we either choose it to be diffusive, 
i.e. = — KaVaM, or we take it to be proportional to £ x A, where A is the vector 
potential of the mean fleld .B = V x A. T he latter expression follows from the recent 
realization ( Hubbard fc Brandenburg 2011 ) that terms involving £ should not occur in 
the expression for the flux of the total magnetic helicity. This will be referred to as the 
alternate quenching model. 
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Figure 1. Saturation magnetic energy for different magnetic Reynolds numbers with closed 
boundaries and diffusive fluxes (solid line) and without (dashed line), as well as the alternate 
quenching formalism (dotted line). 



Without diffusive magnetic helicity fluxes (kq. = 0), quenching is not alleviated and 
the equilibrium magnetic energy decreases as Rcj^^ (Fig.[T|). We find that diffusive mag- 
netic helicity fiuxes through the mid-plane can alleviate the catastrophic a quenching 



466 



Simon Candelaresi & Axel Brandenburg 



and allow for magnetic field strengths close to equipartition. The diffusive fluxes ensure 
that magnetic helicity of the small-scale field is moved from one half of the domain to 
the other where it has opposite sign. With the alternate quenching formalism we ob- 
tain larger values than with the usual dynamical a-quenching-even without the diffusive 
flux term. The magnetic energies are however higher than expecte d from simulations 
( Brandenburg fc Subramanian 2005t iHubbard fc Brandenbure 2011 ) , which raises ques- 
tions about the accuracy of the model or its implementation. 



3. Behavior of the aO, dynamo 

In this second part we address the implications arising from adding shear to the system 
and study the symmetry properties of the magnetic field in a full domain. The large scale 
velocity field in equation (|2.1[) is then U = (0, Sz, 0), where S is the shearing amplitude 
and z the spatial coordinate. We normalize the forcing amplitude ao and the shearing 
amplitude S conveniently: 

Ca = J- Cs = TJ, (3.1) 

with the smallest wave vector ki . 

First we perform runs for the upper half of the domain using closed (perfect conductor 
or PC) and open (vertical field or VF) boundaries and impose cither a symmetric or an 
antisymmetric mode for the magnetic field by adjusting the boundary condition at the 
mid-plane. A helical forcing is applied, which increases linearly from the mid-plane. The 
critical values for the forcing and the shear parameter for which dynamo action occurs 
are shown in Fig. [51 




Figure 2. Critical values for the forcing amplitude Ca and the shear amplitude Cs for an 
aSl-dynamo in 1-D mean-field to get excited. The circles denote oscillating solutions, while the 
squares denote stationary solutions. 

Imposing the parity of the magnetic field is however unsatisfactory, since it a priori 
excludes mixed modes. Accordingly we compute the evolution of full domain systems with 
closed boundaries and follow the evolution of the parity of the magnetic field. The parity 
is defined such that it is 1 for a symmetric magnetic field and —1 for an antisymmetric 
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one: 



Es-Ea 



S/A 



r [B{z)±B{-z)f dz, 
Jo 



(3.2) 



Es+ Ea jo 
with the domain height H. In direct numerical simulations Bx{z) and By{z) are horizontal 
averages. The field reaches an antisymmetric solution after some resistive time tj-os = 
1/ {r]kf) (Fig. [3]), which depends on the forcing amplitude Cq. To check whether symmetric 
modes can be stable, a symmetric initial field is imposed. This however evolves into a 
symmetric field too (Fig.|4]), from which we conclude that it is the stable mode. 
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Figure 3. Parity of the magnetic field versus 
time for a random initial field in f-D mean- 
field. 



Figure 4. Parity of the magnetic field ver- 
sus time for a symmetric initial field in f-D 
mean-field. 



The mean- field results are tested in 3-D direct numerical simulations (DNS); Figs. [5] 
and [S] The behavior is similar to the mean-field results. The preferred mode is always 
the antisymmetric one and the time for hipping increases with the forcing amplitude Ca ■ 
This is however very preliminary work and has to be studied in more detail. 
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Figure 5. Parity of the magnetic field versus 
time for a random initial field in 3-D DNS. 



Figure 6. Parity of the magnetic field versus 
time for a symmetric initial field in 3-D DNS. 



4. Conclusions 

The present work has shown that the magnetic helicity fiux divergences within the 
domain are able to alleviate catastrophic quenching. This is also true for the fluxes 
implied by the alternate dynamical quenching model of lHubbard &: Brandenburg] (|201l[ ). 
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However, those results deserve further numerical verification. Further, we have shown 
that, for the model with magnetic helicity fluxes through the mid-plane, the preferred 
mode is indeed dipolar, i.e. of odd parity. Here, both mean-fleld models and DNS are 
found to be in agreement. 
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Discussion 

Sacha Brun: Is there a reason that your system prefers antisymmetric solutions? It 
seems linked to your choice of parameters. 

Simon Candelaresi: So far we do not see a reason for that. But we see a parameter 
dependence of the transition time. We will look at the growth rate of the modes indepen- 
dence of the parameters. This will give us some better clue if also mixed or symmetric 
modes are preferred. 

Gustavo Guerrero: Is there a regime where the advective flux removes all the mean 
fleld out of the domain? 

Simon Candelaresi: If the advective flux is too high the magnetic field gets shed before 
it is enhanced, which kills the dynamo. So, there is a window for the advection strength 
for which it is beneficial for the dynamo. 



